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Abstract. We consider a car-following model described by a delay difference equation 
and give its exact solutions that present propagation of a traffic jam. This model is 
a discrete-time version of the delayed optimal-velocity model; in the continuum limit, 
we recover the delay differential equation for this model and the exact solutions as 
well. We then work in the ultra-discrete limit, obtaining a delay cellular-automaton 
model, which successfully inherits the solutions. Also the dispersion relation for the 
present solutions suggests that a quick response of drivers does not always result in 
fast dissolution of a traffic jam. 
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1. Introduction 

Study of traffic flow relies in no small part on tlieoretical/matliematical modeling 
because of the difficulty in extensive measurements and designed experiments, and 
predictions from simulation and analytical study have particular importance in this 
held p. 

A car-following model is described in general by an equation of the form: 

Xn{t) = G {hn{t - r)) (hnit) != Xn+l{t) - Xn{t)). 

Xn{t) denotes the position of the nth car at time t and a positive constant r is a delay 
in time, which corresponds to the driver’s reaction time. After this latency, each car 
reaches the velocity G{hn(t — r)) determined by the headway — r), i.e., a distance 
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to the next car in front; hence G'(-) is nowadays referred to as the optimal-velocity (OV) 
function. 

In the earlier studies of traffic ffow, G. F. Newell proposed the car-following model 

with 

G(h) = Vo (^1 - exp - L) j J , 

where Vq, 7 and L are positive constants, and found an exact solution for r = 0 |2]. In 
this case, the equation can be transformed into a linear equation and thus allows for 
superposition of the shocks each of which present a traffic jam propagating upstream at 
different speeds. 

Subsequently, G. B. Whitham obtained an elliptic solution for Newell’s model in the 
generic case r 7 ^ 0 |3]. This presents nonlinear waves beyond Newell’s /mear solutions, 
including a solitary-wave in the case that the modulus of the elliptic function goes to 1 . 
He also suggested that the solitary-wave solution be stable in the range of parameters 
for stable linear waves. Both Newell’s model and the Kac-van Moerbeke system reduce 
to the same equation if one assumes a steady-prohle solution with the phase velocity 
l/2r; the Kac-van Moerbeke system is known as one of the integrable lattice systems 
|1]. This fact motivated Whitham to derive a solitary-wave solution of Newell’s model 
after a long calculation. (It is notable that these results consist with the theory of 
nonlinear waves [5].) 

Bando et al. proposed an alternative car-following model described by a second- 
order differential equation Xnit) = a[K(h„(f)) — Xnit)] with the sensitivity a (a positive 
constant) and the OV function, 

I/(/i) = tanh(h — c)-I- tanhc, (1) 

where c > 0 is the inflection point of the tank function. They referred to the model as 
the optimal-velocity (OV) model [6l [7]. (Note that replacing Xn{t -\- r) with -|- aXn 
for small r = 1/a, one formally recovers the OV model.) Now V(h) is the most popular 
choice for the OV function and is actually suitable for mathematical treatment. In |H1|9], 
they showed exact solutions for the car-following model with the OV function ([ 1 ]), 

An(t) = V(hn(t - r)), 

which they call the delayed optimal-velocity (OV) model. 

The car-following model is totally asymmetric in the lattice variable, and which 
is one of the crucial properties of traffic ffow models. However, in the above works, a 
steady-prohle solution is chosen so that the delay in time should compensate for the 
asymmetric lattice shift. Gonsequently, all the solutions obtained, given by elliptic 
functions, have the same phase velocity of l/ 2 r independently of the form of the OV 
function. An alternative way to exactly solve a traffic-how model is to assume a shock- 
wave solution. 
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Inspired by Hirota’s method m to find the soliton solutions of soliton equations, 
in m they reduced the delayed OV model expressed in 

hn{t + t) = V{hn+l{t - t)) - V(hn(t - t )) 

= tanh(h„+i(t — r) — c) — tanh(h„(t — r) — c), 
to two coupled equations, and thus found a shock-wave solution. It is distinct from the 
previous results introduced above that their solution includes a free parameter which can 
control the shock velocity. They also considered the same scenario for Newell’s model, 
finding a shock solution with velocity 1/r; this solution also includes a free parameter 
but does not allow for the variation of the velocity. 

In this work, we solve exactly the discrete delayed OV model firstly introduced in 
[12] and thereby integrate the previous results of m and [l2]. In Section 2, we derive 
the discrete delayed OV model and give exact results. We also consider the continuum 
limit as the discrete delayed OV model reduces to the delayed OV model. In Section 
3, we consider the ultra-discrete limit of the results obtained in the previous section. 
In this limit, the model is transformed into a cellular automaton; thus we obtain a 
delay cellular automaton with an exact solution. Section 4 is devoted to summary and 
conclusion. 

2. Discretization of the delayed OV model and its exact solution 

Changing the variable as 

gnit) = tanh - c), (3) 

in (l2|), we transform the delayed OV model into an algebraic differential equation; 

gn{t) = (l - gn{tf) {gn+i{t - t) - gn{t - t)) , (4) 

where the transcendental function tanh disappears. Then we discretize (jl]) and obtain 
its exact solution assuming the open boundary condition that an infinite number of cars 
run on the road, i.e., n G Z. 

2.1. Discretization of the delayed OV model 

First we shall review the relationship between (jl]) and the modified Lotka-Volterra (mLV) 
equation given by 

rjit) = Tjit) (1 - arjit)) {rj+i{t) - rj_i(t)), (5) 

where j G Z and a G M. 

If one assumes a traveling-wave solution to (jlD, i.e., gnif) = G(0) with 0 = t + 2nT, 
(jl]) yields 

G'w = (1 - G{ 4 >r) {G(4> + t)-G(<I>-t)), (6) 

where ' stands for the derivative. Meanwhile, changing the variable as rj{t) = 
(1 -I- fj(f))/(2a), from (I5|) we have 

4afj(f) = (l - {fj+iit) - fj-iit)). 


( 7 ) 
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Then, we assume the traveling-wave solution, Vj = with := t/(4a) + jr, to ([7]) 
and reach the same equation as ([6|): 

= (1 - R{^f) {R{^P + r) - R{^ - t)) . 


We thus see that both (jl]) and ([5]) reduce to the same equation ([6]) if a traveling-wave 
solution is assumed. 

Here we introduce the discrete modified Korteweg-de Vries (mKdV) equation [T8] . 


+ ..t 1 + ^-1 


av 


t+i 


^ 1 — avj 


( 8 ) 


where 5 G M. This is known as a discretization of the mLV equation (|5]). (We use 
superscript t as the discrete time variable.) Actually, if we let nj = rj{—6t) then (jS]) 
reduces to ([5]) in the continuum limit where 6 regarded as the discrete time unit tends 
to 0. 

In the same manner as above, we change the variable as Vj = {l+ h*)/(2a) in (jS]) 
and have 


1-27 

-I 

7 




]) = (1 - - (1 - (9) 


where 7 := 5/(4a). Assume the traveling-wave solution, nj = H(4/) with T := f -|- mj 
(m G Z>o), and ([9]) yields 

+ 1) - v{^)) = (1 - H(^)) (1 + + 1)) + 1 + m) 

- (1 - ■(/(T + 1 )) (1 + V{^)) I/(T - m). 

Here we introduce an alternative variable and then consider = U{^) with 
^ := t + 2mn. If U{^) satishes flTU]) . we have 

([/($ + 1) _ um) = (1 - t/(<I))) (1 + U(<!> + 1)) 17($ + 1 + m) 

7 (11) 

- (1 - [/($ + 1 )) (1 + [/(d))) [/($ - m). 


Since ul^ = U{^) ($ = f-|-2mn), we may replace U{^ + l+m) with and U{^ — m) 

with consequently obtain the discrete delayed optimal-velocity (OV) model: 

- <) = (1 - 0 (1 + <+T' - (1 - (1 + <) (12) 

Note that flT^ possesses the asymmetry in the lattice variable n which is required for 
the car-following model. 

Let us see that the discrete delayed OV model flTTD is the discrete version of the 
delayed OV model ([ 1 ]). We restrict ourselves to 7 > 0, for 7 corresponds to the discrete 
time unit. If there exists a smooth function ^„(f) that satishes gnilt) = we have 

- -^(^n(t + 7) -^n(t)) = (l-^n(t))(l + ^n(t + 7))^n+l(t-r + 7) 

7 (13) 

- (1 - Quit + 7 )) (1 + gnit)) gnit - r). 
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where we let r := 7717 . Consequently, taking the limit 7 +0, we successfully recover 

dl]) from ([13]). 

Accordingly, the discrete variable corresponding to the headway hn(t) should be 
dehned from (|3]) as 


It 1 1 1 + “n 

K = c+-^og- 




(14) 


2.2. Exact solution of the discrete delayed OV model 

In order to hnd shock solutions of flT^ . we assume 
, M + NL^K^ 


(15) 


" l + L‘/<'“ 

where A" 7 ^ 1, L 7 ^ 1, M and N {M 7 ^ N) are constants to be determined. Since 0151) 
can transform into , we only consider the case of L > 1. Substituting 

flT^ into flT^ . we obtain 

{LK - 1)M^ - 2L™(L - 1)M + L^{L - 1)A - {LK - 1) = 0, 

L^{LK - l)N‘^ - 2{L - 1)KN + {L - 1)KA - L^{LK - 1) = 0, 

(L + 1)(LA: - 1)MN - {LK + 1){L- 1)(M + N) 

+ {LK + 1){L- 1)A - (L + 1){LK - 1) = 0, 
where A := (1 — 27 ) 77 . Elimination of M and N from ffT 6 |) yields the two cases: 

K = L-2, (17) 

L-1- 47(L"*+i - 1) 


(16) 


K = 


(18) 


L(L-1-47(L-L-™))' 

These are the dispersion relations for the traveling-wave solutions 0151) . 

Remember that corresponds to ([3|) and that the distance h„(f) of cars must be 
positive, and accordingly we shall hnd the solution in the range of 

— tanhc < < 1. (19) 

In the case of (ITTD . from flT6|) we obtain two shock solutions; however these he out of 
the range flT^ . 

In the case of OlSp . together with M and N determined by 016p we obtain the two 
shock solutions: 

(1 - 47 ) (L - 1 ) 1 + K^K-m-i 


ul = l- 


u„ — —1 + 


27(1 - L-"^) 
L-1 


1 -h K^L^ 

1 + k^k-^ 


( 20 ) 


( 21 ) 


27 (L - L-™) 1 -h K^L^ 

Figure [H shows the graphs of an exact solution for the discrete delayed OV model: the 
headway h\^ given by fflTl) and ([2T|). This illustrates a shock wave propagating upstream; 
the shock front corresponds to the tail of a traffic jam. Alternatively, (1201) will show an 
upside-down shock propagating upstream, which illustrates dissolution of a traffic jam; 
the shock front corresponds to the head of a traffic jam. 
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Figure 1. Exact shock solution for the discrete delayed OV model: the headway hl^ 
da of nth car where is given by m with parameters c = 1 , L = 1 . 1 , 7 = 0 . 2 , 
m = 3 at time t = —100, t = 0, t = 100. 


The parameter values used in the simulation are chosen to satisfy the conditions 
discussed below. In general, the range of a shock wave is bounded by the limit values 
at inhnity. Hence we have only to investigate these values in order to hnd the range of 
the solution. 

At hrst, we consider Since L > 1 and m G Z>o, the upper limit of fl2U]) is 

smaller than 1 if and only if 1 — dy > 0 is satished. Now we assume 1 — dy > 0. Since 
L > 1, m E Z>o and 0 < y < 1 /d, taking the limits as n —>■ ±00, we hnd that 
satishes 

(1-47)(L-1) ^ , (1-47)(L 

27(1 - L-™) ^ 27(1 - L-™) 

for all t and n. We introduce the function of L, 


f{L) := 1 


(l-dy)(L-l) 


2 y(l - L-^) ’ 
and then have the derivative 

(1 — dy)(L™'+^ — (m + 1 )L + m) 


f\L) = 


2 yL™+i(l - 


( 22 ) 


This leads to f\L) < 0 for all L > 1. Consequently, if and only if — tanhc < /(I), we 
hnd an L that satishes — tanhc < f{L) and a shock solution that satishes flT^ . From 
/(I) = 1 — (1 — dy)/( 2 my), we explicitly show the condition to have an exact shock 
solution: 


1 

d + 2 m(l + tanhc) 


< y < 


1 

d‘ 


(23) 


(We remark that one can hnd an unphysicalsohition, \u\^\ > 1, unless fl2^ is considered.) 
Assuming ([23]), we can hnd the maximum value Lf oi L such that — tanhc < f{L) is 
satished. Consequently, we have 


/(L) >/(L/) =-tanhc (1 < L < L/). ( 2 d) 

Next we consider the solution dZO- As well as the above discussion, we hnd the 
lower limit of fl2T]) is larger than —1 if and only if 1 — dy > 0 is satished. We also assume 
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1 — 47 > 0 in this case. Taking the limits as n ^ icxo, we hnd that satishes 

L-^(L-l) , L-1 

27 (L - L-™) 27 (L - L-™) 

for all t and n. Now we let 

r- + - 1 )’ 

and thus hnd that (l2Til ranges from g{L) to g{L~^). Then we see the condition to have 
an exact solution, 

1 1 
4(m + 1) ^ ^ ^ 2(m + 1)(1 — tanhc) ’ 
and that there exists the maximum value Lg of L such that 


- tanhc = g{Lg) < g{L) <g{L ^) < I (1 < L < Lg). (25) 

Since f{L) > —1 or g{L~^) < 1 is equivalent to L —1 — 47 (L —L”™) < 0, we hnd K > 1 
for L < max(L/,Lg). 

The phase velocity of the shock solutions fl20l) and fl2ll) is given by 
log L log L 


U = 


log K log 


L-1-47(L'"+1-1) 

L(L-1-47(L-L-"‘)) 


From this formula, we can see that the phase velocity may attain the maximum at an 
arbitrary m, but however it appears complicated. For simplicity, we will discuss this 
point for the continuum version. 


2.3. Continuum limit of the shock solution 


As described above, 5 or 7 takes the role of discrete time unit and tends to 0 when we 
consider the continuum limit. 

The dispersion relation fflSj) implies that L converges to 1 or 1/iF in the continuum 
limit 7 +0. However, it is apparent that 1/iF is never admissible for (1201) . We hence 

let L = 1 + /37 + 0 ( 7 ^), which gives L"” -)■ as 7 -)■ + 0 . 

Also, let K = 6°, mj = r and t be replaced with t/j to consider the continuum 
limit. (Note that the latter t denotes the continuous time variable.) Thus we have the 
exact solutions of (jlj). 


9n(t) 


± 



/3 


X _|_ ^an+0{t-T) 


2(1 - e-h^) 1 + e 




( 26 ) 


and the dispersion relation, 

/3 - 4(e^^ 


e“ = 


1 ) 


/3 - 4(1 - e-/5-) ’ 

which is common to both solutions. The solution (126|) and the dispersion relation (|27)) 
obtained here coincide with those given in [TT]. Also, from fl 2 ^ we have 

1 


r > 


2(1 + tanh c) 


( 28 ) 
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Note that c > 0 ensures r > 1/4 and consequently that the right hand side of (1271) is 
positive. For small r, the phase velocity U = jS/a may be approximated by a power 
series in r: 


U = 


4^2 


1 

7 + S + 


This suggests that a quick response of drivers to the traffic situation does not always 
yield fast propagation of a traffic jam. We hnd the maximum velocity 


U = 


21og(gf) 


at r = ilog|^. 



Figure 2. Exact shock solution for the delayed OV model: the headway h„(t) (|2^ of 
nth car where u„(t) is given in (1261) with minus sign and parameters c = 1, d = 0-2, 
T = 0.6 at time t = —100, t = 0, t = 100. 


Figure 12] shows an exact shock solution: the headway 

which corresponds to the solution with minus sign in (126|) . As expected, it illustrates a 
traffic jam propagating upstream. 

3. Ultra-discrete limit of the discrete delayed OV model 

In this section, we consider the ultra-discrete limit of the discrete delayed optimal- 
velocity model (1T2|) . In this limit, the dependent variable is discretized and thus we 
obtain a delay cellular automaton from the delay discrete equation. 

3.1. The ultra-discrete limit 

The ultra-discrete limit was hrstly developed to directly relate the soliton equations to 
cellular automata |15] . 

We dehne the ultra-discrete limit or ultra-discretization by the transformation from 
a positive real variable n to a ultra-discrete variable U: 

U = lim elogu, 

£ —^H“0 
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where £ is the ultra-discrete parameter. (In usual, a ultra-discrete variable is denoted 
by the capital letter corresponding to that of the variable.) By linking u with U by 

u = (30) 

multiplication/division and addition in the original variable are transformed respectively 
into addition and ‘max’ in the ultra-discrete variable: 

lim £ log(M ■ v) = U + V, lim e: log(M -f n) = max(17, V). 

£->■- 1-0 £->- 1-0 

Thus, an equation including multiplication and addition but no subtraction transforms 
into the ultra-discrete equation described with addition and ‘max’ in the ultra-discrete 
variable. Since addition and ‘max’ become integer functions being restricted to integers, 
we can regard the ultra-discrete equation as a cellular automaton. 


3.2. Ultra-discrete limit of the discrete mLV equation 

As mentioned above, the ultra-discretization runs only for equations described with 
multiplication and addition. We hence change the dependent variable as hj = n)/(l — 
avj) to eliminate subtraction from ([H]): 

l + (a + 2f'jti ^ + + (31) 

^ 1 aVjXl 1 + 

Then, we introduce the ultra-discrete variable Vj and parameters D and A as h* = 

6 = a = e~^l^., and thus obtain from fl3lD the ultra-discrete mKdV equation, 

1//+' + max(0, - D, -A)- max(0, - A) 

= 1// + max(0,1/^1 - D, -A)- max(0,1/^^ - A), 
where D and A take integer values as well as Vj m. 


3.3. Ultra-discrete limit of the discrete delayed OV model 


We remark that the variable changes, Vj = (l-|-n*)/(2a), n* = n*/(l —anj) and n* = 
done thus far are collected into 




av, — 1 

= tanh 


aVj + 1 


2e 


Taking gnit) = tanh(h„(f) — c) into account, we introduce the ultra-discrete variable H\ 
corresponding to the car distance h\^ such that 

'Hi-C 


Un =tanh 


2e 


where C > 0 is an ultra-discrete parameter corresponding to c, i.e., = H^/{2e) 

and c = C/{f2e). Then, in order to ultra-discretize (IT^ . we change the variable as 
u\^ = {1 + — u\f) which means = (u* — U)/{vfj + 1) and 


( 32 ) 
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Accordingly, from fll2p we have 


(1 - 47)h*-+r' + 1 (1 - + 1 


u: 


t+i 


u. 


t—m+l I -r 

n+l “T 


= U, 


+ 1 


(33) 


The condition fl2^ for there to be an exact shock solution suggests that we should 
introduce the ultra-discrete parameter G corresponding to 7 as 

1 - 47 = (34) 

where G > 0; this implies that we focus on around 7 = 1/4. Consequently, substitution 
of fl3^ and fl34)) into (l33D . we obtain the ultra-discrete delayed optimal-velocity (OV) 
model, 

+ max(0, ' - C - G) - max(0, - G) 

= Hi + max(0, - G - G) - max(0, i//"™ - G), 
in the ultra-discrete limit. 


(36) 


3 . 4 . Exact solution of the ultra-discrete delayed OV model 

By taking the ultra-discrete limit, we obtain exact solutions of the ultra-discrete delayed 
OV model fl5S]l from the solutions fl2U]) and fl?T]l of the discrete delayed OV model. 

According to the above procedure to obtain the ultra-discrete delayed OV model, 
we transform (120|) and fl2T|) into 

L-l- 47(L - L-^) + [L - 1 - 47(L™+^ - 


u„ = 


(1 - 47)(L - 1)(1 + 


(36) 


1 — 47 


and 


«n = 


1 K^Lt-m-1 ’ 

(l-i)(i + a:’^l*-”^) 


L-l- 4:^{L - L-™) -L[L-1- 47(L™+i - 1)]K^L^-"^ 

1 1 + 

~ - 1 1 -h ■ 

' L —1 

(We exploit the dispersion relation f[T5]) .i By use of the function f{L) dehned in 
we have 


(37) 


47 


L-L- 


- 1 = (1 - 47 ) 


l + /(^) 


i-fiLY 

is positive and hence the prefactors in 


(38) 


L-l 

From the condition ([23]), we see that 
and fl57|) are positive. 

Now we introduce the ultra-discrete parameters as L = and K = Since 
L > 1, we note that Q > 0. The dispersion relation flTSl) is transformed into 

1 + L + ■ ■ • + L'"-! 


KL -^{1 - 47 ) + 
1 — 47 


1 -F L H-h T™ 

1 + L + ■ ■ ■ + 


KL 


1 + L + ■ ■ ■ + 















11 


and yields 

max(P — mQ — G, max(0, (m — 1)Q) — max(0, mQ)) 

= max(—G — Q,P — mQ + max(0, (m — 1)Q) — max(0, mQ)), 

in the ultra-discrete limit. Note that G > 0, and one may reduce this equation to the 
ultra-discrete dispersion relation: 

max(Q — G, mQ — P) = 0. (39) 


Then let 

T — T 

- 1 = e®/', 

and this, together with flT8|) . transforms into 

KLe^/^ + (1 - 47)(1 + L + ■ ■ ■ + L'") = L + + ■ ■ ■ + L™. 

Thus we have 

max(P Q B, —G max(0, mQ)) = max(Q, mQ) 

in the ultra-discrete limit. Since G > 0, Q > 0 and m G Z>o, we compare both sides 
and consequently obtain 


B = —P {m — 1)Q. 

Finally we obtain exact solutions of the ultra-discrete delayed OV model (l35il by ultra- 
discretizing (136|) and fimi : 

P* = G-|-G — P-|-(m — 1)Q max(0, {n + 1)P + {t — m)Q) 

— max(0, nP {t — m — l)Q), 

P* = G -|- P — (m — 1)Q -|- max(0, nP + {t — m)Q) 

— max(0, {n -|- 1)P -|- (f — m -|- 1)Q). 

These are surely the solutions of fl3^ as far as the dispersion relation fl39|l holds. 

The following conditions are added so that the headway P* should take positive 
values. The range of L given in (|2T)) or fl25|) determines the range of Q for fITOl) or (ITTH 
respectively in the ultra-discrete limit. From fIMll we have 

1 — 47 


(40) 

(41) 




47W-1' 

Then, let c = G/2£ and we have, together with the above results, the condition, 
G + G-P+(m-l)Q>0, 

for firUj) in the ultra-discrete limit. As well, from fl2^ we have the condition 


G > mQ, 


for (ITTll . 

Figure [3] shows an exact shock solution fITTD . It provides a quite simple expression 
of traffic jam propagation. 
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Hk H^n 





Figure 3. Exact shock solution for the ultra-discrete delayed OV model: the headway 
Hn (IdTI) of nth car with parameters C = 4, P = 3, Q = 1, m = 3 at time t = —20, 
t = 0, t = 20. 

4. Summary and conclusion 

In this work, we consider the discrete delayed OV model described by a delay difference 
equation, where the time variable is fully discretized and the delay is hence given by 
an integer. In the continuum limit, the model reduces to the delayed OV model, one 
of classical traffic-flow models, which is described by a delay differential equation. We 
moreover work in the ultra-discrete limit where the dependent variable is also discretized 
to take integer values, proposing a novel traffic-flow model described by a delay cellular 
automaton, which we call the ultra-discrete delayed OV model. 

We find exact solutions presenting a traffic jam of a shock wave form for the 
discrete delay model considered. Working in the continuum and ultra-discrete limit, 
we also obtain the corresponding solutions to the delay differential equation and the 
delay cellular automaton. To our knowledge, this is the first case of a delay cellular 
automaton with an exact solution. 

The significance of the exact solutions presently obtained lies in that, as well as the 
shock formation of vehicular traffic, we have the velocity of a traffic jam propagating 
in an exact form from the dispersion relation. The shock velocity implies that a quick 
response of drivers to the change of car distance does not always result in fast dissolution 
of a traffic jam. We think that moderate response to the deceleration of cars in front will 
be effective to suppress disturbance of a homogeneous flow. We note that in previous 
works [m [m EH], in contrast with the present study, they give exact results for the 
fundamental diagram, i.e., the density-flow plot to investigate a phase transition from 
homogeneous to inhomogeneous density flow; but however these results do not show the 
detailed form of a traffic jam. 

The exact solutions can also prove, as discussed in [3], that nonlinear waves such 
as a shock and a solitary wave may appear in the range of parameters where the linear 
solution, which presents a uniform flow, exists stably. This means that we can explore 
nonlinear waves beyond the linear stability theory. A detailed discussion on this point 
will be found in our subsequent publications. 

We finally add that the present work is motivated by the fact that the delayed OV 
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model shares the same solitary-wave solution with the Kac-van Moerbeke system, one 
of integrable lattice systems, and we have expected that the soliton theory including 
discretization and ultra-discretization of integrable systems applies well to the present 
model. 
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